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Abstract The generalized Chaplygin gas (GCG) modeH 
in spatially flat universe is investigated. The cosmolog- 
ical consequences led by GCG model including the evo- 
lution of EoS parameter, deceleration parameter and di- 
mensionless Hubble parameter are calculated. We show 
that the GCG model behaves as a general quintessence 
model. The GCG model can also represent the pres- 
sureless CDM model at the early time and cosmologi- 
cal constant model at the late time. The dependency 
of transition from decelerated expansion to accelerated 
expansion on the parameters of model is investigated. 
The statefinder parameters r and s in this model are 
derived and the evolutionary trajectories in s — r plane 
are plotted. Finally, based on current observational 
data, we plot the evolutionary trajectories in s — r and 
q ^ r planes for best fit values of the parameters of 
GCG model. It has been shown that although, there 
are similarities between GCG model and other forms of 
chaplygin gas in statefinder plane, but the distance of 
this model from the ACDM fixed point in s — r diagram 
is shorter compare with standard chaplygin gas model. 

Keywords Dark energy. Generalized Chaplygin gas 
model, Statefinder diagnosis 



1 Introduction 

Since 1998, the type la supernova (SNe la) observa- 
tions have shown that the universe has undergone to 
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the accelerating expans ion phase (IRiess et al.. I Il999 ; 
IPerlmutter et al.~ 1999 ). This fact has also been sup- 
ported by many additional observations, including the 
anisotropy measurements of Cosmic Microeave Back- 
ground (CMB) fr om Wilkinson Microwave A nisotropy 



Probe (WMAP) (jSpergel et al. Il2003l 120071 ). the data 



of Large Scale Structure of universe (LSS) from Slo an 



Digital Sky Survey fSD SS) dTeemark et al.. I l2004allbh 



and X-ray experiments ([Allen et al.. 112004) . The cur 



rent accelerating expansion of universe indicates that 
in addition to the existence of dark matter, which is 
required to explain the galactic dyna mics and the for- 



mation of structures ( Bosma. 1 119811 ) , the universe is 



dominated by an exotic energy component with neg- 
ative pressure, dubbed the dark energy. In an an- 
other word, in the framework of standard cosmology, 
the dark energy (DE) scenario is a theoretical solu- 
tion to explain the accelerating expansion of the uni- 
verse. The combined analysis of cosmological observa- 
tions suggest that the universe consists of about 70% 
dark energy, 30% dust matter (cold dark matter plus 
baryons), and negligible radiation. The cosmological 
constant, whose equation of state is independent of 
cosmic time, is a simple solution of DE problem. How- 
ever, it suffers from two well known problems, namely, 
the fine-tuning and the cosmic coiii cidence problems 



( Copeland. Sami fc Tsuiikawa. 20061 ). In addition to 



cosmological constant, many kinds of dynamical DE 
models, whose equation of state is no longer a constant 
but slightly evolves with time, have been suggested to 
interpret the cosm ic acceleration. The quintessenc e 
( Wetterich. Ill988l). phantom field JCaldwell. I liooj) . 
quintom (lElizalde. Noiiri fc Odinstov. 1120041), Chaply- 



gin g as models (IKamenshchik. Moschella fc Pasauier . 
200ll ). K-esse nce JChiba. Okabe fc Yamaguchi. Il2000l). 



ihyo 

19991) and agegraphic (|Cai. Il2007[) DE models are the 



tachy on field ( Sen, ho lograp hic (jCohen. Kaplan fc Nelson7[ 



examples of dynamical DE model. It is emphasized that 
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the predictions of cosm ological constant model is still statefinder pair {r,s} is defined as 
fitted to the observation ( Jassal. Bagla fc Padmanabhan. ^ 
20041) . Therefore, a suggested dynamical DE model 
should not be faraway from cosmological constant. Be- 
sides the DE models, modifie d gravity theori e s suc h 



3(g-i/2)' 



(1) 



as scalar tensor c osmology (Boisseau et al.. 200ol). 



braneworld models ( Dvali. Gabadadze fc Porrati. Il2008h | 
have been suggested to solve the accelerated expansion 
of universe. 

The Chaplygin gas is one of the candidate of DE mod- 
els to explain the accelerated expansion of the universe. 
The striking features of Chaplygin gas DE is that it 
can be assumed as a possible unification of dark matter 
and DE. The Chaplygin gas plays a dual role at dif- 
ferent epoch of the history of the universe: it can be 
as a dust-like matter in the early time, and as a cos- 
mological constant at the late time. This model from 
the field theory points of view h as been investigated in 



It is clear that the statefinder is a geometrical di- 
agnostic, because it depends only on the scale fac- 
tor. The role of statefinder pair is to distinguish 
the behaviors of cosmological evolution of dark en- 
ergy models with the same values of Hq and qq at 
the present time. Up to now, the statefinder diag- 
nostic tool has been used to study the various dark 
energy models. The various DE models have differ- 
ent evolutionary trajectories in {r,s} plane, therefore 
the statefinder is a good tool to discriminate DE mod- 
els. For example, the well-known ACDM model corre- 
sponds to a fixed point {r = l,s = 0} in {r,s} plane 
Also, the quintessen ce DE model 



Sahni et al.. 


20031) 


Sahni et al.. 


2003; 



Alam et al.~l l2003bl). the inter 



emerges as an effective fiuid associated with D-branes 
( Bordemann &: Hoppe. Ill993) and can also be obtained 



Zhang. 



from the Born-Infeld action (jBento. Bertolami fc Sen 
20041) . The simplest form of Chaplygin gas model called 



( Zhang 



( |Bihc, Tupper fc Viollier, 11200^). The Chaplygin gas acting quintess ence models dZimdahl fc Pavon. I [2004 : 



2005ar) . the holographic da rk energy models 



2005bl : IZhang et al.. 1 120071 ). the holographic 



standard Chaplygin gas (SCG) which has been used to 



dark energy model in non-flat universe ( Setare. J. Zhang fc X. Zhang. | 



expla in the accelerated expansion of universe (jCori. 
20041) . Although the SCG model can interpret the ac- 



2007D, t he phantom inodel (jChang et al.. 1120071) . the 
tachyon ( Shao fc Gui. 112007 ). the agegraphic DE model 



celerated expansion of universe, but it can not explain 
the astrophysical problems such as the structure for- 
mation and cosm ological perturbation power spectrum 



with and without interaction in fiat and non-fiat uni- 

verse (IWei fc Cai. ll2007HMalekiani fc Khodam-Mohammadi. j 



2010[) and the interacting new a gegraphic DE model 



( Sandvik. |[2003) . Subsequently, the SCG is extended 



in flat and non-flat universe (IZhang et al 



2010: 



into the generalized Chaplygin gas (GCG) which could 
construct viable cosmological models. Same as SCG 
model, the GCG model can obtain th e accelerated ex- 
pansion of the universe ( Setare. 2007t ). 

The quantities H = d/a and q = —d/aH^, namely, 
the Hubble parameter and the deceleration parameter, 
are the geometrical parameters to describe the expan- 
sion history of universe, where a is the scale factor and 
dot denotes the derivative with respect to time. It is 
obvious that a > 0{H > 0) indicates the expansion of 
universe and a > 0{q < 0) means that the universe 
is undergoing an accelerated expansion. The various 
DE models give the same value of go at present time, 
therefore, the bubble parameter H (first time deriva- 
tive of scale factor) and the deceleration parameter q 
(second time derivative of scale factor) can not discrim- 
inate the various DE models. For this aim, we need the 
higher order of time derivative of scale fac tor. Using the 



Khodam-Mohammadi fc Malekiani et al.. I 120101 ). are 
analyzed through the statefinder diagnostic tool. The 

SCG model is analyzed i n terms of statefinder with or 

without dust component ( Gorini. Kamenshchik fc MoschellaT^ 
2003). In 2003, the generalized cosmic chaplygin 
gas (GCCG) model is introduced by Gonzalez-Diaz 



( Gonzalez-Diaz. 1 120031) . The interesting features of 
GCCG model is that it can be stable and free from 
un physical behaviors even when the v acuum fluid sat 



isfies the phantom energy condition (jGonzalez-Diaz 
20031 ). Chakraborty, et al. (2007) have performed the 



third order time deriv ative, Sahni, et al. ([Sahni et al 



20031) and Alam et al. l^Alam et al.. Il2003al ) introduced 



the stateflnder pair {r,s} in order to remove the de- 
generacy of Ho and qo of different DE models. The 



statefinder analysis for GCCG model and for particu- 
lar choice of interaction parameter, they have shown 
the role of st atefinder parameters for t he evolution of 
the universe (jChakrabortv et al.. II2007I ). Zhang, et al., 
(2006) proposed a new model to describe the unifica- 
tion of dark matter and DE, namely the new general- 
ized chaplygin gas (NGCG) model and calculated the 
cosmological consequences and statefinder analysis for 
this model. The advantage of NGCG model is that it 
can represent the other forms of dark energy models 
such as quintessence-like and phan tom- like dark energy 



( X. Zhang. We fc J. Zhang. 1120061) 



In this paper, first we study the cosmological con- 
sequences of GCG model, and then examine it by 
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means of statcfinder diagnostic tool. Here, we also 
study the dependency of the cosmological quantities 
and statefinder diagnostic on the parameters of GCG 
model. The paper is organized as follows: In section 2, 
we introduce the GCG model and derive the statefinder 
parameters {r,s} for this model. In section 3, the nu- 
merical results are presented. We conclude in section 
4. 



2 GCG model 

The equation of state of GCG is given by 

P= (2) 

where ^ > and a > are the par ameters of the 



model (jBento. Bertolami fc Sen. 1 120021 ). In the case 
of a = 1, the GCG model is reduced to SCG model. 
Also, for a = it can be reduced to standard ACDM 
model. In the framework of Friedmann- Robertson- 
Walker (FRW) cosmology, using Eq. ([2]) and the con- 
servation equation d{pa^) ~ ~pd{a^), the energy den- 
sity of GCG is written as 



PGCG = Poggg[As + (1 - ^)a-^(i+")]TT^, 



(3) 



where a is the scale factor, Ag — Aj p]^^Q and Poggg 
is the present value of energy density. Using Eq. ([2]) 
and ([3]), the equation of state (EoS) parameter of GCG 
model can be obtained as 



v4,a3(l+a) 



WGGG 



1- A, +^a3(i+") 



(4) 



From Eq. it is clear to show that at the early time 
(a — > 0), the EoS parameter tends to zero {wgcg — >■ 0) 
and at the late time (a cxd): wqcg ~^ ~1j which are 
equal to the EoS parameter of matter and cosmological 
constant, respectively. Therefore, the EoS parameter of 
GCG model is constrained to the interval — 1 < wa < 
and can be viewed as a general quintessence DE model. 
Moreover, Eq. (|4]) shows that at the early time the 
GCG fluid can be interpreted as a CDM {wgcg 
0), and at the late time it mimics the ACDM model 
{wgcg ^ ^1)- Wc can also see that for As > \ the 
EoS parameter wgcg < ^1 and the GCG model can 
cross the phantom divide. 

Since the dark matter and dark energy are unified by 
GCG model, therefore it can be de-composited into two 
components of dark matter and dark energy, as pgcg = 
Pde + Pdm- Also, by assuming the pressureless CDM, we 
have Pgcg — Pde- Considering the evolving density of 
CDM as 



Pdm — POdmO. 



(5) 



it is obvious that the energy density of DE in GCG 
model can be derived as 



Pde 



Pgcg ~ Pdn 



(6) 



Pogcg[As + (1 - ^)a-3(i+")]T^ - 



POdma 



-3 



By assuming that the universe is filled by GCG compo- 
nent (DE+CDM) and baryonic matter component, the 
total energy density is pt = pcCG + Pb- In the case of 
flat universe, the friedmann equation for GCG model is 
written as 



"3~ 



Pt 



[pogcg[As + [l - As)a-^('+"^]Th 



PobO. 



(7) 
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where pob is the present density of baryonic matter. 
Substituting the dimensionless parameters 



POGCG 



8^G 



GCG 



POb 



8ttG 

^GCG + f^b = 1, 



(8) 



in Eq. ([7]), the Hubble parameter is expressed as 

=HlE^{a)^ (9) 
Hl({l - ^t,)[As + (1 - ^)a-3(i+")]T^ + r!ba-3 

where E{a), the the normalized Hubble parameter, is 
defined as 

E{a) = ((1 - nb)[As + (1 - ^)a-3(i+")]TT^ + Q^a'^^ '^'| 

(10) 



Recently, Xu and Lu ( Xu &: Lu. 2010l ). by applying 
the Markov Chain Monte Carlo approach on the latest 
observational data, have constrained the GCG model. 
The observational data that have been used are: the 



constitution dataset (jHicken. 1120091 ) including 397 type 
superno va la (SNIa) , the observational Hubble data 



(OHDjJSimon 



tion (Allen et al 



20051). t he cluster X-ray gas mass frac- 
20081 ). the measurement results of 



baryon acoustic osc illation (BAO) fr om Sloan Digital 
Sky Survey (SDSS) (|Eisenstein. IboOSi ) and T wo Degree 
Field Galaxy Redshift Survey (2dFGRS) |Percival. ' 



20091 ) ■ and the cosmic micro wave background (CMB) 
data from five-year WMAP |Komatsu. I l2009l ). They 



obtained that in the flat universe, the best fit values of 
the GCG model parameters {Ag, a) and the cosmolog- 
ical parameters {Vli,h? , Hq) with their confidence level 



are: Ag = 0.76 



+0.029 
'0.039 



{la) 



+0.034 
-0.046 



{2a), a = 0.033: 



-0.066 
-0.071 
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(la) t"oZ (2a), ^},h^ = 0.0233^°;™?^ (la) ^O;^^^ (2a) 
and Ho = 69.97l2:?8 (1^^) ^3:08 (2^), with minimum 
chi-square Xmm = 519.342. At following, we derive the 
deceleration parameter q and the statefinder pair {r,s} 
for GCG model. 

The deceleration parameter q, which denotes the ex- 
pansion phase of the universe, is given by 



1 = -§^-' 



Re- witting q in terms of E, we have 
I dE 

q{x) = - 1 

E ama 



(11) 



(12) 



Substituting E from Eq. ([TU| in p^ . we obtain the 
deceleration parameter q for GCG model as 



(1 - A,)aa-3(i+") 
"^, + (1-^3)0-3(1+") 



(18) 



For a = or = 1, we obtain {r = 1, s = 0} which 
refers to the statefinder pair of spatially flat ACDM 
model. Departure of a given dark energy model from 
the fixed point {r = l,s = 0} is a criterion for cval- 
uating of this mode l from spatially flat ACDM model 
( Sahni et al.. II2OO3I ) . The importance of the statefinder 
diagnostic is that the current values of the parameters 
s and r can be extracted from the observational data 
of SNAP (Super Nova Acceleration Probe) type exper- 
iments. Therefore, the statefinder diagnostic combined 
with the future SNAP observation can help us to dis- 
criminate between different dark energy models. 



3(l-r!fc)fA, + (l-^,)a"3(i+«) 



c(l -A,)a-3(i+") +3rjf,a- 
2(1 - Qb) (As + (1 - As)a''^'+"'^ ^ 



-3 



1. 



(13) 



Eq. ([T^ explicitly shows the dependence of decelera- 
tion parameter q on the GCG model parameter As and 
a. Now, we derive the statefinder pair {r,s} for GCG 
model. Using the definition of statfinder parameters in 
Eq. (dJ, we have 



H 



3q 



(14) 



Similar to q, the parameter r can be re- written in terms 
of E as 



1 cPE I , dE -i dE 

r(3^)' + -T^ + l, 



Edilnaf E^^dlna E dlna 
also the parameter s can be obtained as 

1 d'^E I 1 ( dE \2 I 3 dE 
_ E d{lnay ~^ E^^dlna) '^Edlna 



3 dE 



(15) 



(16) 



E dlna 



Substituting E from Eq. ^ in (|T5j) and ([161), we ob- 
tain r and s for GCG model as 



9(l-r!6)(l- As) X 



(17) 



(As + (1 - As)a-'^'+"'>) A,aa-3(i+")" / 



3 Numerical results 

In this section we discuss the cosmological consequences 
led by GCG model. For this aim, the evolution of EoS 
parameter of GCG model, wgcg, the deceleration pa- 
rameter q and cosmological evolution of dimensionless 
bubble parameter ,E, are studied. Then we study the 
GCG model by means of statefinder diagnostic point of 



3.1 EoS parameter 

Solving Eq. Q , the evolution of wgcg as a function of 
scale factor for different model parameters As and a is 
shown in Fig. (1). In upper panel, by fixing As = 0.76 
based on observational constrain, we vary the param- 
eter a as 1,0.1,0.01. For a < 1, increasing the pa- 
rameter a leads to a larger value of wqcg- While, at 
a > 1, Wgcg is smaller for larger value of a. Here, we 
see the dual role of GCG model at different epoch of 
the history of universe: It can be assumed as a pres- 
sureless dark matter {w = 0.0) at the early time and 
a cosmological constant with ui = — 1 at the late time. 
In lower panel, by fixing a — 0.033 based on observa- 
tional constrain, we plot wgcg for different illustrative 
values of Ag. For As = 0.0, we see that wgcg = 0.0, 
which refers to the pressureless CDM model. In the 
case of As = 1.0, we have wgcg = ^li which denotes 
the cosmological constant. For < A, < 1, we can see 
— 1 < Wgcg < 0, denoting the time varying equation 
of state of GCG model and representing the general 
quintessence behavior of this model. Also, increasing 
the parameter As leads to smaller value of wgcg ■ 



2(1 - Qb) (As + (1 - ^)a-3(i+")') + 2n 



Ibtt 
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3.2 Deceleration parameter 

Here we calculate the evolution of deceleration param- 
eter, q for a universe dominated by GCG model and 
investigate the dependency of q on the parameters of 
model. In Fig. (2), by solving Eq. p^ . the evolution of q 
as a function of scale factor for different illustrative val- 
ues of model parameters a and is calculated. Here 
we adopt the observational value of fl^h^ as 0.0233. 
First we fix the parameter Ag as 0.76 and vary the pa- 
rameter a as 1.0, 0.1 and 0.01 (upper panel). Then, by 
fixing a = 0.033, we choose the illustrative values for 
As as 0.0, 0.1, 0.2, 0.3 and 1.0 (lower panel). In upper 
panel, one can see that by increasing the parameter a 
q becomes larger for a < 1 and smaller for a > 1. The 
transition from decelerated expansion (q > 0) to accel- 
erated expansion (g < 0) takes place earlier for smaller 
value of a. From this figure we see that q ~ 1/2 at 
the early time (CDM dominated universe) and q tends 
to —1 at the late time (ACDM dominated universe). 
This fact can also be seen in the lower panel for any 
value of As in the interval < As < 1. For As = 0.0, 
we have the constant deceleration parameter q = 1/2, 
corresponding to the CDM dominated universe. For 
the values in the interval < As < 1, q starts from 
1/2 at the early time and leads to —1 at the late time. 
Furthermore, q becomes larger for lower values of Ag. 
Transition from decelerated expansion (g > 0) to ac- 
celerated expansion {q < 0) occurs sooner for higher 
values of Ag. 

3.3 Hubble parameter 

At following, we study the Hubble parameter which 
evaluates the expansion rate of the universe for GCG 
cosmology. Using Eq. (fTO|) . we plot the cosmological 
evolution of E{a) in Fig. (3). First, we fix the coef- 
ficient As = 0.76 and vary the model parameter a as 
1.0, 0.1 and 0.01 (see upper panel). In this case, we 
see that the larger value the parameter a is taken, the 
bigger value the Hubble expansion rate E{a) gets. In 
lower panel, by fixing a, we vary the coefficient As as 
0.0, 0.1, 0.2, 0.3 and 1.0. The case of As = 1.0 repre- 
sents the standard ACDM model and Ag = 0.0 denotes 
the CDM model. Here, one can sec that for larger value 
of As, the Hubble expansion rate E{a) becomes smaller 
at a < 1 and larger at a > 1. Therefore, from the above 
analysis, we find that both parameters As and a can 
impact the cosmic expansion history in GCG model. 

3.4 Statefinder parameters 

The statefinder pair {r,s} for GCG model is given by 
Eqs. PT|) and (HH]), respectively. In statefinder plane. 



the horizontal axis is defined by the parameter s and 
vertical axis by the parameter r. In this diagram, the 
standard ACDM model corresponds to a fixed point 
{r = 1, s = 0}. At the early time, a — >■ 0, the statefinder 
pair {r,s} defined in Eqs. ([TT]) and (fT8|) are reduced as 
follows: 

9 (1 - nt){l - ^)-^A,aa3("+i) 

r = J hi, s = -a, 

2 {1 - - As) — 

(19) 

therefore we see that at the early time, a — > 0, the 
statefinder pair {r,s} for GCG model are: {s = —a, r = 
1}. From Eqs. pT|) and (|18p . we can also obtain 
the statefinder pair {r,s} at the late time, a — > oo, 
as {r = l,s = 0.0}. Hence, the GCG model mimics 
the cosmological constant at the late time. In Fig. (4), 
we show the evolutionary trajectories of GCG model in 
statefinder plane. In upper panel, by fixing As = 0.76, 
we choose the illustrative values 0.1, 0.2 and 0.3 for 
a. While the universe expands, the trajectories of the 
statefinder start from the fixed point {s = — a, r = 1} at 
the early time. The parameter r starts to increase and 
then decares, while the parameter s increases from the 
initial value s = —a at the early time to s = 0.0 at the 
late time. Here, we can easily see that the statefinder 
trajectory is dependent on the parameter a of GCG 
model. Different values of a give the different evolu- 
tionary trajectories in {s, r} plane. The colore points on 
the curves represent the today's values of statefinder pa- 
rameters (sq, tq) and the star symbol indicates the stan- 
dard ACDM model. The distance to ACDM fixed point 
becomes shorter for smaller value of a. We can also see 
that for smaller values of a, the parameter s increases 
and the parameter r decreases. In lower panel, by fixing 
a = 0.033, we plot the evolutionary trajectories in s — r 
diagram for different values of the parameter As- From 
Eqs. Eqs. (fTT]) and (fTSj) . we have {r = l,s — —a} for 
As = 0.0. Therefore, for =0.0 we have a fixed point 
{s = —a, 7' = 1} in s — r plane. Also form Eqs. (|17p 
and (fT8|l one can see that {s = 0, ?- = 1} for As ~ 1.0. 
Hence, in the case of Ag = 1.0, the statefinder param- 
eters {s,r} are coincide to the ACDM fixed point in 
s — r plane. For different values of Ag in the interval 
< < 1, we have the different evolutionary trajec- 
tories in s — r plane (see the right panel). Therefore, 
the parameter Ag affects the evolutionary trajectories 
in s — r plane. The distance to ACDM fixed point 
{s ~ 0.0, r = 1.0} becomes shorter for larger values of 
Ag. The color points represent the today's values of 
the parameters (sq, tq) for different values of Ag. Here, 
we see that the parameter sq increases for larger values 
of Ag and the parameter tq is largest for r = 0.5. In 
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Fig. (5), we plot the evolutionary trajectory in s — r 
plane (upper panel) and q — r plane (lower panel) for the 
best fit observational values: As = 0.76 and a = 0.033. 
In upper panel the evolutionary trajectory starts from 
(s ~ — 0.033, r — 1) at the past time, reaches to the 
(s = — 0.033, r = 1) at the present time (circle point) 
and ended at (s = 0, r = 1) at the future. This behavior 
of GCG model in statefide r plane is similar to NGCG 



model at the early time (IX. Zhang. We fc J. Zhang. 
20061) . where they found that the universe starts from 
the init ial value (r = 1, s = —a) in s — r plane. G orini, 



et al, (jGorini. Kamenshchik &: Moschella. 1 120031) cal 



culated the trajectory of SCG in s — r plane and 
showed that the universe in SCG model starts form 
(s = —1,7' = 1) reaches to {sq ~ — 0.3, ?■ = 1.9) at the 
present time and finally mimics the ACDM model at 
the late time. Therefore the distance of (so7 ''o)in GCG 
model constrained by the above observational value 
from the standard fixed point ACDM (s = 0,r = 1) 
is shorter compare with SCG model. As a similarity, 
we see that for both model the universe mimics the 
ACDM model at the late time. Moreover, the behavior 
of the trajectory is s — r plane is similar for both model, 
where by expanding the universe, r increases to a max- 
imum value then decreases to r — 1 at the late time 
and the parameter s increases forever. In lower panel 
the evolutionary trajectory in q — r plane starts from 
{q — 0.5, r = 1) at the past (note that (g = 0.5, r — 1) 
corresponds to the CDM dominated universe), reaches 
to {q = —0.575, r = 1.26) at the present time and ended 
at {q = —0.965, r = 1) at the future. 



we exhibit the cosmological evolution of E{a) (normal- 
ized Hubble parameter, E{a) = H{a)/Ha). For GCG 
model, both the parameters and a affect the cosmo- 
logical evolution. Finally, we performed the statefinder 
diagnostic tool on the GCG model. Since many cos- 
mological models have been proposed to interpret the 
accelerated expansion of universe, the statefinder diag- 
nostic tool with the parameters r and s which are con- 
structed by higher order derivative of the scale factor 
is needed to discriminate between them. Moreover, the 
present values of r and s can be viewed as a discrimina- 
tor for testing a given dark energy model if it can be ex- 
tracted from observational data in a model-independent 
way. Here we derived the statefinder parameters r and 
s for GCG model and studied the evolutionary trajecto- 
ries of this model in s — r plane. The dependence of the 
evolutionary trajectories and the today's value of {s,r} 
on the model parameters A^ and a has been investi- 
gated. The lower value of a and higher value of As re- 
sult the shorter distance from standard ACDM model in 
s — 7' diagram. Eventually, we plotted the evolutionary 
trajectory of GCG model in s — r and q~r plane based 
on current observational data and found that the dis- 
tance of GCG model from the standard ACDM model 
in s — r plane is shorter compare with SCG model. 
However, both SCG and GCG models have a similar 
trajectories in s — r diagram. Furthermore, the behavior 
of GCG model in statefide r plane is similar to NGCG 



model at the early time (jX. Zhang. We &: J. Zhang. 



I2OO6I) . where the universe expands from the initial value 
(r = 1, s = —a) in s — r plane. 



4 Conclusion 

Summarizing this work, we investigated the generalized 
Chaplygin gas (GCG) model in spatially flat universe. 
Here we studied the cosmological consequences of GCG 
model by calculating the evolution of EoS parameter 
WGCGi deceleration parameter q and cosmological evo- 
lution E{a) of GCG model. In the GCG cosmology, 
the universe starts from the CDM-dominated phase at 
the early time to the DE-dominated universe at the late 
time. The GCG fluid as a general quintessence dark en- 
ergy model can be viewed as a pressureless matter fluid 
(wgcg = 0.0) at the early time and as a cosmological 
constant {wgcg = ^1) at the late time. We also ob- 
tained the deceleration parameter q in GCG model and 
studied the evolutionary treatment of q function 
of scale factor in this model. In GCG model, the pa- 
rameter q starts from the initial value 1/2 at the early 
time (CDM-dominated universe) and converges to — 1 
at the late time (A-dominated universe). Furthermore, 
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x=ln (a) 




In (a) 




Fig. 1 The evolution of EoS parameter of GCG model, 
wgcg, versus 2; = Ina for different illustrative values of 
parameters As and a. In upper panel, fixing As by best fit 
value: As = 0.76, we vary a as 1,0.1,0.01 corresponding 
to black solid line, blue dashed line and red dotted-dashed 
line, respectively. The case of a — 1.0 exhibits the standard 
Chaplygin gas (SCG) model. In lower panel, fixing a by best 
fit value: a = 0.033, As is varied as (black thick solid line), 
0.1 (blue solid line), 0.2 (red dashed line), 0.3 (green dotted- 
dashed line) and 1.0 (black dashed line). The cases of As = 
0.0 and As = 1.0 exhibit the EoS parameter of pressureless 
matter and cosmological constant fiuids, respectively. 



Fig. 2 The evolution of deceleration parameter q in GCG 
model versus a; = Ina for different illustrative values of 
model parameters As and a. In upper panel, by fixing 
As = 0.76, we vary a as 0.2, 0.1, 0.01 corresponding to black 
solid line, blue dashed line and red dotted-dashed line, re- 
spectively. The case of a = 1.0 represents the SCG model. 
In lower panel, by fixing a = 0.033, As is varied as (black 
thick solid line), 0.1 (blue solid line), 0.2 (red dashed line), 
0.3 (green dotted-dashed line) and 1. (black thick dashed 
line). The cases of = 0.0 and As = 1.0 represent the evo- 
lution of q in CDM-dominated and A-dominated universe, 
respectively. 
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15 




In (a) 



Fig. 3 Cosmological evolution of normalized Hubble pa- 
rameter as a function of logarithmic scale factor x = Ina 
in GCG model. In uper panel, we choose the observational 
best fit values: As = 0.76 and Qbh^ = 0.0233 and vary 
the parameter a as 1.0, 0.1 and 0.01 corresponding to black 
solid, blue dashed and red dotted-dashed lines, respectively. 
The case of a = 1.0 exhibit the SCG model. In lower panel, 
by fixing a = 0.033, is varied as 0.0, 0.1, 0.2, 0.3 and 1.0 
corresponding to black thick solid line, blue solid line, red 
dashed line, green dotted-dashed line, black thick dashed 
line, respectively. 
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Fig. 4 An illustrative example for the statefinder diag- 
nostic of GCG model. In upper panel, the evolutionary 
trajectories in s — r plane are plotted, by fixing As — 0.76 
and varying a as 0.1, 0.2 and 0.3 corresponding to black 
solid line, blue dashed line and red dotted dashed line, re- 
spectively. The circle point on the curves show the today's 
value of statefinder parameters (so,ro). The star symbol 
indicates the location of standard ACDM model in s — r 
plane:{s = 0, r = 1}. In lower panel, the evolutionary tra- 
jectories are plotted for different illustrative values of As, 
by fixing a = 0.1. The case of As = 0.0 is related to square 
symbol located at {s = -0.1, r = 1.0}. Note that As = 0.0 
represents the CDM model. The case of yl^ = 1.0 is exhib- 
ited by star symbol at: {s = 0, r = 1.0} which is related to 
ACDM model. The evolutionary trajectories of illustrative 
cases As = 0.25, As — 0.50 and As = 0.75 have been shown 
by black solid line, blue dashed line and red dotted-dashed 
line, respectively. Circle point on the curves denotes the 
today's value (so,ro) in s — r plane. 
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1.035 
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Fig. 5 The statefinder diagrams r{s) (upper panel) and 
r{q) (lower panel) for GCG model. The evolutionary 
trajectories are plotted in the light of best fit result of 
SNe-HOHD+BAO-hCMB, a = 0.033 and As = 0.76. The 
circlepoints on the curves show the today's value (so,''o), 
upper panel, and {qo,ro), lower panel. For comparison, the 
standard ACDM model has been shown by star symbol in 
these diagrams. 
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